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ABSTRACT 

We present deep radio observations of the most distant complete quasar sample drawn 
from the Sloan Digital Sky Survey. Combining our new data with those from literature 
we obtain a sample which is ~ 100 per cent complete down to Si.4ghz = 60/xJy over 
the redshift range 3.8 < z < 5. The fraction of radio detections is relatively high (~ 43 
per cent), similar to what observed locally in bright optical surveys. Even though the 
combined radio and optical properties of quasars remain overall unchanged from z ^ 5 
to the local Universe, there is some evidence for a slight over-abundance of radio-loud 
sources at the highest redshifts when compared with the lower-z regime. 

Exploiting the deep radio VLA observations we present the first attempt to di- 
rectly derive the radio luminosity function of bright quasars at z > 4. The unique 
depth ~ both in radio and optical - allows us to thoroughly explore the population 
of optically bright FR II quasars up to z ~ 5 and opens a window on the behaviour 
of the brightest FR I sources. A close investigation of the space density of radio loud 
quasars also suggests a differential evolution, with the more luminous sources showing 
a less pronounced cut-off at high z when compared with the less luminous ones. 
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1 INTRODUCTION 

During the past years, studies of the properties of Active 
Galactic Nuclei (AGNs) and in particular of their cosmo- 
. logical evolution have become of major relevance within the 
' more general field of galaxy evolution. In fact, it has been 
found that the properties of the central black hole (BH) are 
tightly related to those of the host galaxy (e.g. Magorrian 
et al. 1998; Ferrarese & Merritt 2000; Tremaine et al. 2002; 
McLure & Dunlop 2004), so that the energetic feedback that 
AGN activity can release is now a fundamental ingredient in 
many theoretical models of galaxy formation (Silk & Rees 
1998; Fabian 1999; Granato et al. 2001, 2004; Cavaliere & 
Vittorini 2002; Di Matteo et al. 2005; Cirasuolo et al. 2005b). 

Previous studies in the optical and X-rays - bands 
which are thought to trace the accretion processes onto the 
central BH - have established quasar and powerful AGN 
activity to peak at 2; ~ 2, with a rapid decline at lower red- 
shifts (e.g. Boyle et al. 2000; Ueda et al. 2003; Groom et 
al. 2004); on the other hand, the less powerful sources have 
their major shining phase at z < 1 (e.g. Ueda et al. 2003; 
Hasinger et al. 2005). The behaviour of the AGN evolution 
at higher redshifts has been very uncertain for a long time 
as the relevant observations were biased by selection effects 
and only considered very small numbers of objects. The ad- 



vent of the recent Sloan Digital Sky Survey (SDSS; Fan et 
al. 1999, York et al. 2000) has allowed to properly explore 
the high redshift Universe up to 2 ^ 6 (Fan et al. 2001a,b; 
2004). These studies have confirmed the presence of a cut 
off in the space density of quasars at z ~ 2. 

Even though AGNs that show radio emission are only 
a small fraction of the total population (Sramek & Weed- 
man 1980; Condon et al. 1981; Marshall 1987; Miller, Pea- 
cock & Mead 1990; Kellermann et al. 1989), they repre- 
sent an important subsample as the radiation at centimetre 
wavelengths is unaffected by dust obscuration and redden- 
ing. Therefore, studies of the evolution of radio-active AGNs 
provide a less biased view of the behaviour of massive BHs 
and accretion processes onto them as a function of cosmic 
time. 

Several studies have addressed the evolutionary trend 
of radio loud sources from the local Universe up to high z 
(Dunlop & Peacock 1990; Toffolatti et al. 1998; Jackson & 
Wall 1999; De Zotti et al. 2005). Shaver et al. (1996, 1999) 
argued for a drop in the space density of flat spectrum radio 
quasars by more than a factor 10 between z ~ 2.5 and z ~ 6. 
However, a re-analysis of such sources (Jarvis & Rawlings 
2000) suggests a more gradual (factor ~ 4) decline, decline 
which is backed up by the work of Jarvis et al. (2001) on 
steep spectrum sources in the same redshift interval. A lu- 
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minosity dependent cut-off, with a decrease in space density 
less dramatic for the most luminous radio sources, has been 
claimed by Dunlop (1998) and confirmed by other recent 
studies (i.e. Vigotti et al. 2003; Cirasuolo et al. 2005a). 

Unfortunately, the process (cs) responsible for the for- 
mation of radio jets that mark the class of radio loud ob- 
jects are still poorly understood. The mass of the central 
BH could play an important role in shaping the transition 
between the population of radio loud (RL) and radio quiet 
(RQ) AGNs. As recently pointed out by many studies, RL 
sources seem to have the BHs confined to the upper end of 
the BH mass function, whereas the BHs in RQ quasars ap- 
pears to span the full range in BH mass (Laor 2000; McLure 
& Dunlop 2002; Dunlop et al. 2003; Marziani et al. 2003; 
McLure & Jarvis 2004; Metcalf & Maghocchetti 2006). Fur- 
thermore, the analysis of a large sample of local low luminos- 
ity AGNs drawn from SDSS showed the fraction of galaxies 
hosting a radio-loud AGN to be a strong function of BH 
and stellar mass (Best et al. 2005). However, the point is 
still controversial, and some authors claim no evidence for 
any relation between radio power and mass of the central 
BH (Oshlack et al. 2002; Ho 2002; Woo & Urry 2002a,b; 
but see the dissenting view of .Jarvis & Mclurc 2002 and 
McLure & Jarvis 2004, who ascribe the lack of correlation 
reported by these latter authors as due to selection effects 
such as Doppler Beaming and orientation). 

In the light of the above discussion, the present work is 
aimed at exploring the radio properties of the highest red- 
shift quasars. The main goal is to investigate if the physical 
conditions in the eaxly stages of galaxy formation can favour 
or prevent the formation of relativistic jets and also to test if 
the radio loudness in quasars exhibits some dependence on 
cosmic epoch. For this purpose, wc performed deep radio ob- 
servations of a sample of high redshift quasars selected from 
SDSS. The optical sample is presented in Section 2, while 
radio observations and the radio properties of the sample 
are respectively described in Sections 3 and 4. By exploiting 
this unique sample we derive the radio luminosity function 
in Section 5 and investigate the behaviour of the space den- 
sity of high redshift QSOs as a function of redshift in Section 
6. Our discussion and conclusions are presented in Section 7. 
Throughout this paper wc adopt the "concordance" cosmol- 
ogy, consistent with the Wilkinson Microwave Anisotropy 
Probe data (Bennett et al. 2003), i.e.: = 0.3, Oa = 0.7 
and Ho = 70 km s~^. 



2 THE SAMPLE 

The optical quasar sample we consider in this work was 
drawn from the Sloan Digital Sky Survey (SDSS) and rep- 
resents the largest, complete, high-redshift quasar sample to 
date. 

This dataset (Fan et al. 2001a) consists of 39 bright 
{i* < 20), high-redshift (3.6 < z < 5) quasars observed in 
the Fall Equatorial Stripe at high Galactic latitude, covering 
an area of 182 sq. deg. High redshift quasar candidates were 
selected using colour cuts by exploiting the five photometric 
bands {u* , g* , r* , i* , z*) available from SDSS observations. 

The quasar continuum is assumed to be a power law 
with a slope Uo, i.e. fu oc . The continuum magnitude, 
ABi4Bo, is defined as the rest-frame AB magnitude at A = 



1450 A corrected for Galactic extinction. For a power law 
continuum, the AB1450 magnitude can be converted into the 
rest-frame Kron-Cousins B band magnitude: 

B = AB1450 + 2.5ao log(4400/1450) + 0.12, (1) 

where 4400 A is the effective wavelength of the B band filter, 
and the factor 0.12 is the zeropoint difference between the 
AB magnitude system and the Vega-based system (Schmidt 
et al. 1995). 

The mean optical spectral index for the sources in the 
sample is Oo = 0.79 (Fan et al. 2001a), even though Penter- 
icci et al. (2003) - making use of near-IR observations (J, H 
and K band) found a slightly flatter slope, Oo ^ 0.5. 

Due to the optical limit i* < 20 adopted to select these 
sources, the final sample is only representative of bright 
quasars with absolute magnitudes in the range —28.5 < 
Mb < -26. 



3 VLA OBSERVATIONS 

Out of the 39 quasars which form the optical sample, 5 
of them have radio counterparts in the FIRST (Faint Im- 
ages of the Radio Sky at 20 cm) survey (Becker, White & 
Helfand 1995) with relatively high 1.4 GHz fluxes: 1.4 < 
>S'i.4GHz (mjy) < 10. These sources have an offset between 
radio and optical positions of less than 0.4 arcsec, except for 
the faintest one (J2322-0052) which has an offset of 0.8 arc- 
sec, still compatible with the positional accuracy of 1 arcsec 
of the FIRST survey for fluxes ~ 1 mJy. Their structure is 
point-like, except for J0210-0018 which is the brightest of 
our objects at 1.4 GHz with 5'i.4ghz ~ 9.75 ± 0.1 mJy; the 
FIRST image of this source reveals a small component with 
a flux of 2 mJy at a distance of 10 arcsec from the centre of 
the optical position, probably indicating the presence of an 
extended structure (with a proper linear size of ~ 65 kpc). 

Other 15 sources of the original Fan et al. (2001a) sam- 
ple were observed by Carilli et al. (2001), as a part of an 
extensive program on high redshift quasaxs at radio wave- 
lengths. These sources were observed at 1.4 GHz for ~ 2 
hours each, reaching a theoretical rms noise (arms) of the or- 
der of 20 — 30 /iJy. Only 4 sources were detected at > 3 fTrms 
and other 3 sources at 2.5 < Crms < 3. Considering the exact 
coincidence of the radio positions with the optical ones, in 
the following analysis we will take as true detections even 
these latter three sources with > 2.5 (Trms. 

In order to constrain the radio properties of the whole 
optically selected sample, we then performed deep 1.4 GHz 
observations with the Very Large Array (VLA) of the re- 
maining 19 sources. These observations were made in the 
BnA configuration in October 2003 and January 2004; the 
correlator setup was such that the total bandwidth was 100 
MHz with two orthogonal polarizations. Every source was 
observed for about 2 hours, divided into short scans span- 
ning a large range in hour angle, in order to have a reason- 
ably good coverage of the u-v plane. The data reduction and 
analysis were performed within AIPS; standard phase and 
amplitude calibration were applied, and subsequently these 
were refined with self-calibration using background sources 
in the field of view. The dirty maps were cleaned and then 
restored using Gaussian beams of about 5 arcsec. The the- 
oretical noise level has been roughly achieved in most of 
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Figure 1. Images at 1.4 GHz of the 4 high redshift quasars detected at > 3<Trms- The FWHM of the Gaussian restoring beams are shown 
in the corners of all frames. 



the sources, except for J2306+0108 (<Trms = 65 /iJy) due to 
side-lobe confusion by a bright source in the field. 

Four quEisars (see Figure 1) have been detected at 
> SfTrms- Two of these sources, J0244-0108 and J0059+0003, 
are relatively bright (S'i.4GHz = 612 ± 20 /ttJy and 
S'i.4GHz = 461 ± 29 /iJy, respectively), while the other 
two have fluxes in the range 150 < S'i.4ghz (/^Jy) ^ 180 (see 
Table 1). Another source (J0204-0112) has been detected at 
~ 2.5 fTrms with a flux S'i.4GHz = 53 ± 20 jiSy. The re- 
maining 14 sources arc undetected and their fluxes will be 
treated as upper limits in the following analysis. 



4 REDSHIFT DISTRIBUTION AND 
COMPLETENESS 

Combining our VLA observations with the ones made by 
Carilli et al. (2001) and with the information stemming from 
FIRST we obtain 17 radio detections out of the 39 quasars 
(~ 43 per cent) constituting the whole optical sample. 

The redshift distribution of these radio-detected sources 
is shown in Figure 2 as a solid histogram. For a compari- 
son. Figure 2 also shows the redshift distribution of the 39 
quasars of the optical sample. It is worth noting that the 
fraction of radio detections is relatively high (~ 40 per cent) 
at z ^ 3.8 but it drops to ~ 10 per cent in the redshift range 



3-6 i$ 2 < 3.8. We are not aware of any specific selection ef- 
fect that can justify this rather peculiar trend. One possible 
cause of such finding would be a significant contamination 
of the sample by young radio sources, whose spectrum can 
be inverted up to tens of GHz and whose density is expected 
to increase at higher redshifts. 

Despite the results of a Kolmogorov-Smironv (KS) test, 
which show that the redshift distribution of radio detected 
sources and that of the underlying parent population are not 
statistically different {Pks ~ 0.3), in order to avoid any un- 
foreseen bias in our analysis, in the following we will always 
discuss results obtained for z > 3.8 and only comment on 
the effects of including sources from the 3.6 ^ z < 3.8 bin. 

As shown in Figure 3, the radio follow up of the optical 
sample in the range 3.8 ^ z < 5 is nearly complete down to 
5'i.4GHz ~ 60 juJy. The only exception is a single upper limit 
for the source J2309-0031 from CariUi et al. (2001), which 
has a 2 a detection at Si.4ghz = 88 ± 44 /i,Jy with relatively 
high noise. Only one detected source (J0204-0112) lies under 
the flux limit of 5'i.4ghz ~ 60 /iJy and it is our less secure 
detection at a 2.5 a level. Note that pushing the limit down 
to z < 3.8 would have the advantage of including another 
source with 5i.4ghz > 60 fJ-Jy, but would imply dealing with 
three more upper limits rather than true detections. 

The radio flux limit of iSi.4GHz ~ 60 fj,3y al- 
lows to reach completeness down to logioPi.4GHz ~ 23.6 
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Table 1. Properties of the 19 quasars observed at 1.4 GHz with the VLA. Redshift and optical data are from Fan et al. (2001a). Upper 
limits to the undetected sources have been defined as 5i.4ghz = 2 (Trms- 
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Figure 2. Redshift distribution for the radio detected quasars 
(solid histogram) compared with the one for the 39 sources in the 
Fan et al. (2001a) dataset (dotted histogram). 

W Hz~^ sr~^over the entire redshift range considered in this 
work. This threshold in radio power of course depends on the 
adopted radio spectral index*, an- Unfortunately, to our 
knowledge no quasar of the Fan et aJ. (2001a) sample was 
observed at radio frequencies other than 1.4 GHz. We were 
therefore unable to directly measure the radio spectral index 
of the sources in our sample and in the following analysis the 

* Throughout this work the radio flux density is defined as 



value qr = 0.5 has boon adopted. Since in the considered 
redshift range observations at 1.4 GHz select the rest frame 
spectrum where the flat core component dominates (~ 8 
GHz), it is possible that these high-z sources have a rather 
flatter radio spectral index than what assumed. The effects 
of the adoption of a flatter radio spectral index (i.e. q_r = 0) 
on our results will be discussed throughout this paper. 

It is also interesting to look at the combined radio and 
optical properties of these high z quasars. In Figure 4 the 
distribution of radio-to-optical ratios^ (^1.4) for the radio 
detected sources in the high redshift sample is shown as a 
solid histogram while the one derived by using the upper 
limits is shown as a dotted line. For a comparison, we also 
plotted as a dashed lino the RI 4 distribution obtained by 
Cirasuolo et al. (2003b) for quasars with z < 2, renormal- 
izcd to match the number of sources in our sample. Quali- 
tatively, there is an overall agreement between the Rl,4 dis- 
tributions at high and low redshifts, at least for Ji*.4 ^ 10; 
both distributions show a region of stoop transition between 
the RL and RQ regime. This statement is also confirmed by 
the KS test which gives high probability Pks ^ 0.6 that 
the two distributions are not significantly different. In the 
radio loud regime, instead, there is some evidence for an 
excess of RL sources in the high-z dataset. In fact, there 
are 5 sources with S'i.4GHz > 1 mjy in such RL tail: two 
of them have R1.4 ~ 10, close to the RL/RQ transition re- 
gion and the other 3 have Ri,4 > 100. By using the RIaghz 
distribution from Cirasuolo et al. the expected number of 
sources with Ri,4 ^ 10 (and 5'i.4ghz > 1 mJy) in the same 
range of redshift and luminosity is ~ 2, suggesting an over- 
abundance of RL sources at the highest z. However, due to 
the limited statistics (only 5 detected objects), the former 

t The 4 parameter is defined as the ratio between the rest 
frame radio luminosity at 1.4 GHz and the optical luminosity in 
the B band. 
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Figure 3. Redshift versus radio power for the 39 sources in the 
complete sample. Filled circles show detections obtained from this 

work, open circles arc from Carilli ct al. (2001) and filled squares 
are detections from FIRST. Long and short down arrows repre- 
sent upper limits from this work and from Carilli ct al. (2001), 
respectively. The solid line describes the selection effect due to a 
radio flux limit 5i.4ghz = 60 /iJy, while the vertical line shows 
the cut at 2 = 3.8 (see text for details). 



estimate is still compatible with the observed number fig- 
ures. In fact, a KS test applied between the distribution of 
RL sources (i?^.4 > 10) at « > 4 and that obtained at lower 
redshifts gives probabilities 0.2 < Pks ^ 0.5 (depending 
on the adopted value of an) suggesting that the two dis- 
tributions are not statistically different. This result remains 
unchanged even if - in order to increase the statistics - lower 
values of the radio-to-optical ratios {Ri.i > 2— 3) are consid- 
ered. Therefore, even if there is some evidence for an excess 
of radio-loud sources at high z, we cannot make any strong 
statement. Larger samples are definitely needed to clarify 
this issue. 



Figure 4. Distribution of radio-to-optical ratios for radio de- 
tected sources (solid histogram) and undetected ones (dotted his- 
togram), for which upper limits in radio power have been used. 
The dashed line is the ^ distribution. 
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Table 2. Binned radio luminosity function in the range 3.8 < 
z <5 for Q.m = 0.3, Ha = 0.7 and Hq = 70 km s"! Mpc"!, as 
plotted in Figure 5. 



5 LUMINOSITY FUNCTION 

Even though small, the sample of radio detected quasars de- 
scribed in the previous Section is sufficient to perform some 
statistical studies and derive for the first time an estimate 
of the radio luminosity function (RLF) of optically selected 
quasars at high redshift. By considering both the optical 
and radio completeness limits of the sample (i* < 20 and 
S'i.4GHz > 60 /iJy, respectively), we can derive the RLF by 
using the classical 1/Vmax method (Schmidt 1968). For each 
source we evaluated the maximum redshift at which it could 
have been inchided in the sample, Zmax = min(zmax, z°ax), 
where z^^^ and z"^^ respectively correspond to the values 
implied by the radio or the optical limiting flux densities. 
The incompleteness related to the optical multicolour selec- 
tion has been taken into account by applying to each object 
in the sample a detection probability as provided by Fan et 



al. (2001a) and used by Fan et al. (2001b) to compute the 
optical luminosity function. 

Figure 5 shows the RLF obtained via the 1/Knax 
method for sources with iSi.4GHz > 60 /xJy in the redshift 
range 3.8 < z < 5 (see also Table 2). As described in the 
previous Section, within this redshift range the radio follow- 
up is almost complete with the exception of a single upper 
limit. In order to test the reliability of our estimate of the 
RLF, we then re-computed it by assuming this upper limit 
to be a real detection. We find that the above assumption 
only affects the RLF in a negligible way as only the faintest 
bin changes by a factor less than 5 per cent. This gives us 
evidence for the stability of our results. We also computed 
the RLF assuming a flatter radio spectral index [or = 0). 
The main effect is to translate the RLF towards lower radio 
powers by a factor 2.5 without affecting its shape, since 
the determination of the maximum accessible volume (Knax) 
is driven by the optical cut. 
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Figure 5. Radio luminosity function for sources in the range 
3.8 < 2 < 5 with Si.4GHz ^ 60 fijy (solid dots). For a compar- 
ison, as a solid line we show the RLF obtained by De Zotti ct 
al. (2005) for flat spectrum radio quasars, calculated at the mean 
2 of the bin. Short dashed and long dashed curves illustrate the 
RLFs obtained by Dunlop & Peacock (1990) for flat spectrum ra- 
dio sources by respectively assuming a pure luminosity evolution 
and a luminosity dependent evolution model. 

For a comparison, Figure 5 also shows the high z ex- 
trapolation of the RLF for flat spectrum radio quasars as 
recently inferred by De Zotti et al. (2005). The evolution- 
ary model derived by De Zotti et al. for the population of 
radio sources has been obtained by fitting number counts 
and redshift distributions of several radio selected samples 
at different frequencies. An overall good agreement between 
our determination of the RLF and the one derived by De 
Zotti et al. (2005) is found over about two decades in radio 
luminosity. 

In Figure 5 we also compare our RLF estimates with the 
predictions by Dunlop & Peacock (1990) for fiat spectrum 
radio sources. The dashed and long dashed curves show the 
RLFs obtained by respectively assuming a pure luminos- 
ity evolution and a luminosity dependent evolution model, 
translated to the concordance cosmology adopted in this pa- 
per. We note a moderate agreement with observations only 
for logioPi.4GHz < 25 W Hz^^ sr^^in the case of luminosity 
dependent evolution of radio sources. However, it is worth 
noticing that the Dunlop & Peacock RLFs are completely 
unconstrained at these high redshifts, so they have to be 
merely considered as a pure extrapolation. 

It is however important to notice that a direct com- 
parison between the RLF obtained from our high redshift 
quasars and models for the different radio populations as de- 
rived by De Zotti et al. (2005) or Dunlop & Peacock (1990) is 
not straightforward. These models in fact include all the ra- 
dio sources regardless of their optical properties, while our 
RLF is only representative of bright {Mb < —26.5) opti- 
cal quasars. For a meaningful datarto-model comparison we 
should than take into account the (unknown) contribution 



from fainter FRII sources and quasars. The problem is how- 
ever not as bad as it looks since - as Cirasuolo et al. (2003b) 
have shown - even though with a large scatter, optical and 
radio luminosities are broadly correlated so that the con- 
tribution of these "missing" faint sources should typically 
affect the RLF at luminosities lower than those probed by 
our work. 

As we will discuss in more detail in Section 6, we have 
also compared our findings with extrapolations of the RLF 
obtained at lower z by Cirasuolo ct al. (2005a) and Willott et 
al. (2001). These extrapolations are not able to reproduce 
the high redshift data presented in this paper since they 
assume a negative evolution which is too quick for 2; > 2 — 3. 



6 SPACE DENSITY 

The number of quasars in our sample only allowed us to com- 
pute the RLF in a single redshift bin. Moreover, as described 
in the previous Section, comparisons with evolutionary mod- 
els for the different radio populations is biased by the bright 
optical cut in our sample and the lack of information on the 
radio spectral indices. Therefore, the best way to explore 
the cosmological evolution of radio quasars and shed light 
on the behaviour of radio loudness up to the highest acces- 
sible redshifts is to look at the evolution of an integrated 
quantity such as their space density. 

Figure 6 shows the comparison between space den- 
sities of radio quasars as computed at different redshifts 
and for different ranges of optical luminosities. The solid 
square in the Figure represents the space density obtained 
from the high redshift sample (3.8 < z < 5) presented in 
this work by considering sources with Mb < —26.5 and 
lo gloPl.4GHz > 24.4 W Hz"^ sr~^ The optical cut is due to 
the apparent magnitude limit i* < 20 of the optical sample. 
The cut in radio power has instead been chosen in order 
to have an unbiased comparison with the space densities 
computed at lower redshifts. For z < 2.2, the space den- 
sity was derived by using a combined sample from the 2dF 
Quasar Redshift Survey (Croom et al. 2004) and the Large 
Bright Quasar Survey (LBQS; Hcwett et al. 1995). Both 
these surveys have been cross-correlated with the FIRST 
dataset (Cirasuolo et al. 2003a, 2005; Hewett et al. 2001) 
in order to obtain samples of quasars complete in radio 
down to a flux limit of S'i.4GHz = 1 mjy, which translates 
into logioPi.4GHz > 24.4 WHz"^ sr~^at z < 2.2 (open cir- 
cles). In order to minimize selection effects in the compari- 
son of the space densities of high and low redshift quasars, 
we then applied to the low redshift sample the same opti- 
cal {Mb < —26.5) cut of the 2 ~ 4 dataset. It is worth 
noting that some selection effects can arise from the differ- 
ent selection techniques applied to construct the 2dF-LBQS 
samples and that of SDSS quasars. However, at the bright 
magnitudes considered in this paper these effects should be 
negligible. 

The filled dots in Figure 6 show the space density ob- 
tained from the combined 2dF-LBQS-FIRST sample as de- 
scribed above. It is clear from the Figure that the drop in 
the space density of optically bright RL quasars between 
z ^ 2 and z ^ 4.4 is not very pronounced although er- 
rors are large. This result remains mostly unchanged when 
the space density is computed over the entire redshift range 
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Figure 6. Space density of radio quasars with logioPi.4GHz > 
24.4 WHz"'^ sr^'^computed at different redshifts and for differ- 
ent optical luminosities (sec text for details). The filled square 
represents the space density for those sources in the high redshift 
sample (3.8 < z < 5) presented in this paper {Mg < —26.5), 
while the filled dots show the space density for the combined 
2dF-LBQS-FIRST sample at lower redshifts derived by applying 
the same optical and radio cuts as above. The space density of 
optically fainter sources {Mg < —24.4) from the combined 2dF- 
LBQS-FIRST sample as obtained by Cirasuolo et al. (2005) is 
shown at the different redshifts as open circles. 



3.6 < < 5. In fact, no other source with logioPi.4GHz > 
24.4 W Hz^^ sr^^ is present in the range 3.6 < z < 3.8, and 
the effect of sampling a slightly larger volume is minimal. 

To quantify the decline in the space density of optically 
bright radio quasars we apply a linear regression to the three 
data points with Mb < —26.5, obtaining a value for the 
slope of —0.07 ± 0.11, entirely consistent with no evolution. 
If we instead use a flatter radio spectral index {or ~ 0) 
for the high redshift sources, the space density at a ~ 4.4 
drops by a factor ~ 1.5 leading to a slightly steeper slope 
of —0.14 ± 0.13. Therefore, we conclude that, independent 
of the adopted radio spectral index, the drop in the space 
density of optically bright RL quasars between z ~ 2 and 
z ~ 4.4 is at most a factor ~ 1.5 — 2. 

The above trend is in agreement with recent results 
from Vigotti et al. (2003). By using a sample of 13 radio 
quasars, these authors showed the decline in their space den- 
sity to be approximately a factor 2 between z ^ 2 and 4. 
This is close to our findings, even though the objects used 
by Vigotti et al. (2003) arc about one magnitude brighter 



in optical {Me 



< 



-27.5) than those in our sample, while 



the radio selection in the two datasets is roughly the same 
(logioPi.4GHz > 24.6 WHz-i sr-i). 

It is interesting to note that the behaviour of bright 
quasars scorns to differ from that of fainter sources which 
instead present a sharper decline toward higher redshifts. 
Figure 6 shows (open circles) the space density at low red- 
shifts {z < 2.2) of fainter {Mb < —24.4, completeness limit 



of the combined 2dF-LBQS-FIRST dataset for z < 2.2) op- 
tical quasars. As already pointed out by Cirasuolo et al. 
(2005a), this space density shows a decrement of a factor ~ 2 
already in the narrow redshift range z ~ 1.7 and z ~ 2.2. In 
fact, by fitting with a straight line the values of the space 
density at z > 1.7 we find a slope of —0.6 ± 0.3. This trend, 
when compared with the behaviour of the optically brighter 
sources which only show a moderate decline up to ^: ~ 4.4, is 
suggestive of a differential evolution of radio-active sources 
of different optical luminosities. Certainly, the lack of in- 
formation on the space density of fainter sources at higher 
redshifts {z > 2.2) does not allow any definitive conclusion, 
but our results indicate the cosmological evolution of radio 
activity in quasars to be a function of their optical power. 



7 DISCUSSION AND CONCLUSIONS 

We have presented deep radio observations of the most dis- 
tant complete quasar sample drawn from the Sloan Digital 
Sky Survey. Combining our deep VLA observations with the 
ones performed by Carilli et al. (2001) and also with the 5 
detections from FIRST, we obtained ~ 100 per cent com- 
pleteness down to 5i.4GHz = 60 /iJy over the redshift range 
3.8<z< 5. 

The fraction of radio detections is relatively high (~ 43 
per cent), similar to what observed locally for bright optical 
surveys such as the Palomar Bright Quasar Survey (KcUcr- 
mann et al. 1989) and the LBQS (Hewett et al. 2001). A 
comparison between the i?i.4 distribution of these high red- 
shift radio quasars with the one derived at z < 2 by Cira- 
suolo et al. (2003b) suggests that the combined radio and op- 
tical properties of quasars might remain overall unchanged 
from « ~ 5 to the local Universe. However, even though 
the shape of the R*a distribution is roughly preserved over 
cosmic time, there is some marginal evidence for a slight 
over- abundance of radio loud sources at high z when com- 
pared with the low redshift samples (see Figure 4), even 
though not statistically significant due to the small number 
of sources. Furthermore, it is worth noting that the adoption 
of a flatter an, for high z objects would reduce their radio 
power and shift them towards lower values of Rt,4, therefore 
somehow reducing the fraction of purely RL {Ria) sources. 

An interesting hint to shed some light on the above issue 
comes from comparisons with the parent optical population. 
The space density of bright optical quasars {Mb < —26.5) 
at z ~ 2 and 2 ~ 4.4 is po ~ 2 x 10"^Mpc~^ and 
po ~ 1-5 X 10~^Mpc~^, respectively (see Fan et al. 2001b, 
2004). The ratio between the space densities of radio sources 
with logioPl.4GHz > 24.4 W Hz"^ sr"^(as plotted in Fi gure 
6) and the po of the total optical population is therefore 
0.025 ± 0.01 at z ~ 2 and 0.15 ± 0.1 at z ~ 4.4. These fig- 
ures have been obtained for a radio spectral index or = 0.5 
but, as shown in section 6, the results arc not expected to 
exhibit great variations by adopting a flatter slope. Again, 
this suggests that at high redshifts the probability of having 
RL sources is enhanced with respect to that at lower red- 
shifts. However, we stress once more that the statistics we 
dealt with in this work is very poor and further data and 
larger samples are needed in order to have a more robust 
answer. 

Even though a detailed investigation of this phenomenon is 
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outside the possibility of the present data, in a very quali- 
tative way we can relate the suggested excess of RL sources 
at high z as compared with the lower redshift regime with 
changes in the accretion rates with cosmic time. The phys- 
ical conditions of the primordial massive galaxies hosting 
quasars and the availability of a larger amount of gas in these 
early stages could in fact allow super-Eddington accretions 
and eventually favour the formation of powerful radio jets. 

We also attempted the first direct estimate of the radio 
luminosity function of bright quasars at 2: > 4. Exploiting 
the deep radio flux limits obtained through VLA observa- 
tions, we wore able to trace the RLE down to logioPi.4GHz ~ 
23.6 WHz^^ sr^^. It is worth noticing that the transition 
region between the FR I aud PR II population occurs at 
logioPi.4GHz ~ 24.4 WHz^'^sr^^. Therefore, the unique 
depth - both in radio and optical - of this high redshift 
quasar sample allows us to completely explore the popula- 
tion of optically bright ER II quasars up to z ~ 5 and fur- 
thermore opens a window on the behaviour of the brightest 
FR I sources. 

Einally, close investigation of the RL quasar space den- 
sity at different redshifts is suggestive of a differential evo- 
lution for the two populations of optically faint and bright 
objects. The more luminous sources in fact show a less pro- 
nounced cut-off at high z - with a drop in their space density 
of only a factor ~ 2 between ~ 2 and z ~ 4.4 - when com- 
pared with the less luminous ones. Even though the lack 
of information on the behaviour of optically faint quasars 
at « > 2.2 does not allow any definitive conclusion, our re- 
sults indicate the cosmological evolution of radio activity in 
quasars to be a function of their optical power. 
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